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ABSTRACT 


An  analytical  approximation  to  the  image  theory  model  is  developed  for  predicting  the 
acoustical  pressure  field  in  a  wedge-shaped  ocean.  This  thesis  is  a  continuation  of  the  ongoing 
development  of  this  model.  Previously,  limitations  has  restricted  the  source  and  receiver  to  be  placed 
in  only  an  upslope/downslope  configuration.  This  thesis  removes  these  limitations  and  allows  the 
source  and  the  receiver  to  be  placed  in  cross-slope  configurations. 
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INTRODUCTION 


Current  historical  events  have  necessitated  a  change  from 
studying  the  propagation  of  sound  in  deep  water  to  studying 
the  propagation  of  sound  in  shallow  water.  With  the  fall  of 
the  Soviet  Union  and  the  rise  of  Mid-Eastern  countries. 
Antisubmarine  Warfare  (ASW)  is  more  and  more  likely  to  occur 
in  shallow  water.  Shallow  water  provides  a  much  more  complex 
environment  for  ASW  prosecution,  and,  unfortunately,  has  not 
had  as  much  research  as  in  deep  water. 

Image  theory  is  one  method  that  is  being  developed  to 
predict  sound  propagation  in  shallow  water.  This  method  is 
more  accurate  than  ray  tracing  methods  because  it  is  not 
limited  to  high  frequencies.  A  number  of  thesis  conducted  at 
the  Naval  Post  Graduate  School  were  involved  with  applying 
image  theory  to  a  wedge-shaped,  shallow-water  ocean.  The 
first  thesis  provided  the  groundwork  for  the  initial 
development  of  the  image  theory.  The  most  recent 
[Nassopoulos,  1990]  derives  a  simplification  involving  a 
doublet  approximation  and  incorporates  all  of  this  information 
into  one  manageable  computer  program  [Nassopoulos,  1992]. 

This  thesis  assimilates  a  more  complicated  cross-slope 
approximation  into  the  doublet  analysis  which  allows  the 
source  and  receiver  to  be  placed  anywhere  in  the  wedge.  A 
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computer  program  has  also  been  developed  to  show  that  the 
approximations  are  correct  within  the  s^t  of  assumptions. 
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II.  BACKGROUND 


Research  on  the  prediction  of  sound  propagation  in  a 
wedge-shaped  ocean  has  been  directed  in  three  major  areas:  the 
parabolic  equation  approximation,  the  adiabatic  normal  mode 
theory,  and  the  image  theory. 

The  parabolic  equation  can  be  used  as  a  range  dependent 
underwater  acoustic  propagation  model.  It  basically  replaces 
the  Helmholz  equation  with  a  one-way  parabolic  approximation 
which  will  generate  the  acoustic  field  as  an  initial  value 
problem.  The  original  version  of  the  parabolic  equation, 
developed  in  1977  by  Tappert,  contained  restrictions  on  the 
maximum  elevation  angle  at  the  source.  These  limitations  were 
reduced  by  Claerbout  and  Green  and  finally  eliminated  with  the 
use  of  a  higher  order  parabolic  equation  [Collins,  1987]. 
Because  the  parabolic  equation  is  cylindrically  symmetric, 
sloping  bottoms  cause  a  problem.  One  attempt  to  solve  this 
problem  has  been  attempted  by  sectioning  off  the  bottom  into 
a  series  of  range- independent  regions  and  then  applying  the 
normal  parabolic  approximation  to  each  section  [Lee  et  al, 
1983].  Although  this  method  increases  the  accuracy,  it  still 
does  not  alleviate  the  problem  completely.  A  more  recent 
approach  to  this  problem  has  been  attempted  by  using  a  rotated 
parabolic  equation  [Collins,  1990]  which  shows  good  results. 
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Normal  mode  theory  is  a  range-independent  approach.  This 
approach  has  been  used  for  open  ocean  predictions,  but  for  a 
wedge  shaped  ocean,  a  range  dependent  approach  is  needed.  For 
this  type  of  approach  it  was  shown  [Pierce,  1965]  that  an 
approximation  to  the  normal  mode  theory  can  be  utilized  by 
performing  an  adiabatic  separation  of  the  depth  and  range 
coordinates  in  the  wave  equation.  This  approach  has  been 
applied  to  the  wedge-shaped  ocean  [Graves  et  al,  1975]  but  was 
limited  to  only  small  slopes  and  non-penetrable  bottoms.  An 
exact  solution  for  the  ideal  wedge,  where  boundaries  are 
perfectly  reflected,  was  developed.  [Buckingham,  1987]  A 
continuing  problem  for  adiabatic  mode  theory  is  that  it  can 
not  explain  the  transition  from  the  guided  mode  to  the 
evanescent  modes  at  cutoff  [Jensen  et  al,  1980].  One  hybrid 
approach  that  has  been  pursued  to  account  for  this  is  the  use 
of  the  parabolic  equation  [Pierce,  1982]  which  is  solved 
asymptotically  subject  to  adiabatic  mode  initial  conditions. 
Another  method  is  based  on  tracking  of  local  plane-wave 
spectra  [Arnold  et  al,  1984]  which  produce  intrinsic  mode 
fields.  Another  method  uses  spectral  representation  based 
on  the  theory  of  characteristic  Green's  functions  [Kamel  et 
al,  1983]. 

Another  approach  in  solving  the  wedge-shaped  ocean  problem 
is  image  theory.  This  theory  predicts  the  pressure  field  from 
a  series  of  projected  images  [Coppens  et  al,  1980].  A  further 
description  of  this  approach  will  be  provided  in  the  following 
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chapters.  The  basic  geometry  and  mathematical  analysis  has 
been  developed  and  compared  to  the  parabolic  equation  model 
with  favorable  results.  The  upslope/downslope  case  has  been 
approximated  [Nassopoulos ,  1992]  by  use  of  Taylor  Series 
approximations  resulting  in  a  substantial  savings  in  computer 
computational  time.  Although  less  research  time  has  been 
placed  into  the  image  theory  approach,  it  shows  high  potential 
as  a  viable  way  to  predict  acoustic  pressure  fields  in  a  wedge 
shaped  ocean. 
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III.  DEVELOPMENT 


The  basic  goal  of  this  research  is  to  approximate  the 
images  by  a  sum  of  doublets.  Although  a  full  description  of 
the  theories  will  not  be  incorporated  into  this  thesis,  the 
following  can  be  used  as  a  general  overview. 

A.  IMAGE  THEORY 

Image  theory  stems  from  the  idea  that  if  a  source  is 
radiating  near  a  surface,  then  an  image,  out  of  phase  with  the 
source  and  an  equal  distance  away  from  the  surface,  is 
produced.  In  an  isospeed  wedge-shaped  ocean,  numerous  images 
are  produced  [MacPherson  et  al,  1966].  These  images  form  a 
circular  pattern  around  the  point  where  the  surface  meets  the 
bottom  and  interact  with  a  series  of  images  each  possessing  an 
image  of  the  surface  and  bottom  of  the  wedge.  The  images  are 
numbered  sequentially  around  the  upper  half  of  the  circle 
starting  with  1  for  the  source,  2  for  the  image  and  so  on 
around  the  circle.  The  lower  half  of  the  circle  is  labeled  in 
a  similar  fashion.  The  pattern  for  the  sign  convention  along 
with  the  numbering  of  the  images  is  shown  in  Fig.  1. 
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Figure  1 

Each  image  contributes  to  the  total  pressure  field  at  a  point 
in  the  wedge.  For  each  image,  the  distance  to  the  receiver 
and  the  total  reflection  coefficient,  which  is  the  product  of 
all  the  reflection  coefficients  along  the  path  from  the  image 
to  the  receiver,  can  be  determined.  The  total  pressure  is  the 
sum  of  the  pressures  from  all  of  the  images.  A  simplified 
three  dimensional  view  at  Fig.  1  shows  that  the  shoreline  acts 
as  the  perpendicular  axis  to  the  circle  of  images.  See  Fig. 


B.  DOUBLET  RADIATION 

The  images  on  each  side  of  the  pressure  release  surface 
can  be  grouped  together  in  pairs  to  form  acoustic  doublets: 
Images  1  and  2  form  a  doublet,  3  and  4  form  another  doublet, 
and  so  on  around  the  upper  and  lower  half  of  the  circle.  Each 
doublet  consists  of  two  simple  sources,  vibrating  at  the  same 
frequency  but  180  degrees  out  of  phase.  See  Fig.  3. 


Geometrical  Presentation  of  Doublet 
Fieurc  3 

It  can  be  shown  that  in  the  far  field  the  pressure  field 
resulting  from  this  doublet  using  a  far  field  approximation  is 


P(r,0,  t)  =j  —  sin  ( —  kj^dsind)  (1) 

JT  2 

Where  d  is  the  separation  of  the  source  and  image  in  the 
doublet,  r  is  the  distance  from  the  field  point  to  the 
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midpoint  of  the  doublet,  and  B  is  the  angle  formed  by  r  and  a 
perpendicular  bisector  of  the  doublet.  it  should  be  noted 
that,  if  the  polarities  were  switched,  the  pressure  field 
would  change  sign. 

To  incorporate  doublets  into  the  wedge  problem.  Eg.  1 
needs  to  be  modified  to  include  the  reflection  coefficients. 
Nassapolous  developed  the  following  equation  for  pressure 


P(r,e,  t)  =-j  [  (x-iir)  cos(^k^dsind) 
+j  (X'*"'!')  sin(^k^dsind)  ] 


(2) 


Where  x  ^  are  the  cumulative  reflection  coefficients  for 

the  negative  and  positive  images  respectively. 


C.  DOUBLET  ANALYSIS  OF  THE  IMAGE  THEORY 
1.  General  Description: 

The  purpose  of  this  thesis  is  to  extend  Nassapolous' 
doublet  approximation  to  the  image  theory  model  to  incorporate 
cross-slope  propagation.  The  basic  assumptions  remain  the 
same:  the  sound  speeds  and  densities  are  constant  in  both  the 
wedge  and  the  bottom,  the  water-air  boundary  is  a  pressure- 
release  surface,  and  the  slope  is  constant  and  shallow 
[Nassapolous,  1992].  Nassapolous  developed  the  case  for 
upslope  and  downslope  propagation.  Figure  4  shows  how 
doublets  are  grouped  together:  The  third  and  fourth  upper 
images  are  grouped  together  with  the  first  and  second  lower 
images  (n=l) ,  the  fifth  and  sixth  upper  images  are  grouped 
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together  with  the  third  and  fourth  lower  images  (n=2)  and  so 
on  around  the  circle  until  all  the  doublets  are  paired. 


Reprcsentadoa  of  the  Groupinfs  oi  Aeousiicil  Doublets 
Figure  4 


The  source  and  image  for  the  first  surface  reflected  path  form 
the  neutral  doublet  (n=0)  and  is  an  exception  to  this  grouping 
rule.  Similar  to  doublet  radiation,  where  each  doublet 
contributes  to  the  pressure  field,  now,  each  pair  of  doublets 
is  considered  the  basic  contributor  to  the  total  pressure 
field.  The  number  of  paired  doublets  used  in  this  total 
pressure  calculation  can  be  determined  by 


N= 


360 


1. 


)  '2 


(3) 


0  is  the  angle  between  the  surface  and  the  bottom.  To  ensure 
that  the  pressure  field  is  accurate,  a  wedge  angle  should  be 
chosen  such  that  all  of  the  images  are  grouped  into  pairs  of 
doublets.  To  insure  that  there  is  closure  for  the  set  of 
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images,  360/^  mus^  be  a  multiple  of  four.  In  calculating  the 
total  pressure  field,  the  pressure  from  the  neutral  doublet 
will  first  be  developed,  followed  by  the  pressure  field  from 
each  of  the  paired  doublets. 

All  distances  used  are  scaled  according  to  Rscaiea 
~^re«i/^*  Two  bottom  types  are  examined  in  this  research,  a 
fast  bottom  and  a  slow  bottom.  For  a  fast  bottom,  the  scaling 
distance  X,  occurs  where  the  lowest  mode  attains  cut  off.  For 
a  slow  bottom,  the  lowest  mode  does  not  attain  cut  off,  but  a 
scaling  distance  X,  is  used  for  convenience  in  calculations. 
The  definitions  of  these  scaling  distances  are 
Fast  bottom! 


2sin0^tanP 

_ _ n _ 

2tan0ptanp 

cos0^=-^ 

*-2 


(4) 


Slow  bottom; 


2 tan0^ 
n 


tanP 


23in0stanP 
COS0  =-^ 

*  q 


(5) 
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(]  is  the  wedge  angle,  the  subscript  l  refers  to  the  fluid  in 
the  water,  and  the  subscript  2  refers  to  the  bottom 
[Nassopoulos ,  1992]. 

2.  Neutral  Doublet  (n=0) 

The  neutral  doublet,  as  stated  above,  comprises  of  the 
source  and  upper  image.  The  pressure  field  from  the  neutral 
doublet  crosses  no  bottom  boundaries  therefore,  no  reflection 
coefficients  need  to  be  considered.  See  Fig.  5. 


G«Jnieincii  Reprejcnotionof  Nctb»1  Doublet 
Hfurt  5 


d  is  the  distance  between  image  1  (the  source)  and  image  2,  S 
is  the  angle  formed  between  the  surface  and  the  plane 
containing  the  receiver  and  the  shoreline,  and  R2  is  the 
perpendicular  distance  from  the  shoreline  to  the  receiver. 
Taking  Eq.  1,  that  was  developed  for  doublet  radiation  without 
reflection  coefficients,  and  applying  it  to  the  new  cross 
slope  geometry  gives 

Pj.j Msln  (■ik.dsJno,')  (6) 
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where  P is  the  pressure  field  from  the  neutral  doublet,  R'o 
is  the  distance  between  the  receiver  and  the  midpoint  between 
the  neutral  doublet  and  a\  is  the  angle  formed  between  R and 
the  surface. 

To  reduce  Eq.  6,  the  following  substitution  is  made: 

-|=i?iSinY  (7) 


7  is  the  angle  formed  between  the  surface  plane  apd  the  plane 
containing  the  source  and  shoreline,  and  R,  is  the 
perpendicular  distance  from  the  shoreline  to  the  source.  See 
Fig.  6. 


Source 

Geometrical  Representation  of  R|  and  Y 
Figure  6 

Since  S  is  assumed  to  be  a  small  angle,  sin7  can  be 
approximated  by  y  which  reduces  Eq.  6  to 

Po'=J— , sin  (k^F^ysinaJ)  (8) 

A  further  substitution  can  be  made  by  applying  the  geometry  in 
Fig.  7. 
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0«otneirh.al  KcpfcteiHiibwol  andR*^ 

I'igtne  9 

The  side  common  to  both  triangles  shows  that  RjsiniS  = 

R'oSina'o-  Since  S  is  small,  the  small  angle  approximation  can 
be  utilized  to  form  the  equation  sina =  Rj^/R'^. 

Substituting  this  into  Eg.  8  and  assigning  r =  k,R,R2/R'o 
gives  the  final  pressure  for  the  neutral  doublet 

3.  Pressure  Field  From  the  Nth  Pair  of  Acoustic  Doublets 

The  paired  doublets  can  be  split  into  two  groups,  one 
from  the  upper  set  of  doublets  P'„*  and  one  from  the  lower  set 

of  doublets  P'„'-  The  denotes  the  cross  slope  geometry, 

the  subscript  n  denotes  which  paired  doublet  is  being 

calculated  and  the  ”+”  and  the  denotes  the  upper  and  lower 
doublet  respectively.  Figure  8  illustrates  the  geometry  of 
these  upper  and  lower  doublets. 
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'Receiver 

Plane 


nth 

upper 

doublet 


X 

P* - / - ' 

7  Receiver  / 

nth  \  ^ 

lower 

doublet  1 

CieotneiTy  of  Upper  and  Lower  Doublets 
Figure  8 

The  pressure  from  the  upper  doublet  is 


R  '*  2 


j(wt-kRj') 


R 


(10) 


and  the  pressure  from  the  lower  doublet  is 

P/‘  =  -(-l)"j-^sin(-|iCidsino/)  dD 

Rj  2 

The  (-1)  and  the  -(-1)  account  for  the  different 
orientations  of  the  top  and  bottom  doublets. 

The  total  pressure  from  all  of  the  paired  doublets  is 


N 

n-l 

This  pressure  could  be  calculated  by  a  computer  program,  but 
it  would  take  a  relatively  long  time  to  compute.  To  shorten 
the  computational  time,  and  also  research  for  an  analytical 
approximation  a  Taylor  series  expansion  is  used  for  RW 
R'„,  sina'„‘,  and  sina'„‘. 


a.  Derivation  R\*  and  R'^- 

The  geometry  for  an  upper  doublet  image  is 


in  Fig.  9. 


Geometrical  Repiesenudon  of  Upper  Doublet 
From  this  geometry  Rgure9  , 


and 


(R^*)  (2nP+6) 


therefore  by  direct  substitution 


Rn=\lRi^*R2^^Y^^-2R^R^cos{2n^+b) 


By  the  same  reasoning,  for  the  lower  doublet 


shown 


(13) 


(14) 


(15) 
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Rn  ^  V -2i?ii?2COS  ( 2n^  -6 ) 


(16) 


Setting  Yq  zero,  Eq.  15  and  Eq.  16  reduce  to 

Rn=y/Ri"  *^2"  -2R^R^cos  ( 2np  +5 ) 


^n*  ( 2nP  -5 ) 


(18) 


which  match  the  equations  developed  be  Nassapolous. 

(1)  The  Taylor  Series  Expansion  for  R\"  and  R\'- 
Because  R'„*  and  R'„‘  appears  in  the  phase  a  third  order  Taylor 
series  expansion  is  needed.  Talcing  this  expansion  around  5=0 
gives 


Rn-Rno 


6«0 


■  1  52 
2  d62 


6-0 


3! 


db- 


6*0 


(19) 


This  formula  is  the  general  case  for  both  the  upper  and  lower 
doublet.  To  get  the  correct  formula,  R'„  would  be  replaced 
with  R'„*  or  R'„‘  as  needed.  R is  the  distance  to  the  nth 
doublet  if  the  receiver  were  moved  vertically  to  the  surface. 
See  Fig.  10.  Since  the  geometry  is  now  symmetrical 

^no^=-Rno'=i?no=^i?l'-^«2'^^o'-2i?ii?2COS(2nP)  (20) 
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Geotnetrictl  Representation  of  o  and  R^e 
Figure  10  ** 

Taking  the  first  derivative  yields 


dRn  _  i?^i?2sin(2nP+6) 


(21) 


(22) 


d6 


/?^i?2Sin  (2nP) 

R'no 


(23) 


dR'n\ 

db 


-)?ii?2sin  (2/jP) 

rL 


(24) 


Taking  the  second  derivative  yields 
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(25) 


[RnRiR20os{2n^*b)] 
d6^  ~ 

[R^R2Sin(2np+6)  i?^j?^sin(2np-t-5)  ] 
(Rn)^ 


d^Rn 

db^ 


I&-0 


i?j^i?2COS  (2nP) 


(2nP) 


R 


/  3 

JIO 


M»o 


RjR. 


i? 


^[cos(2np)-fiMi21i^l 

'  r>  /  2 


(26) 


(27) 


Through  cancellations  of  the  minus  signs,  the  second 
derivative  for  the  lower  doublet  is  the  same,  so 


d^R'n  ,  _  d^R'n 

db^  db^ 


(28) 


TaKing  the  third  derivative  and  evaluating  at  5=0  yields 


d^R'n 

db^ 


6  *0 


RIR2 


[ 


-3i2j^i?2Cos  {2i2P)  sin  (2np) 

^no 

+  ___  j 


sin(2nP) 

(29) 


d63 


)  6  *0 


R-1^2 


-SRiRjCos  (2nP)  sin(2nP)  .  _ 

[ - - - - —  -  sin(2nP) 

rL 

3R^^R2^sin^  Un^)  ^ 


(30) 


Combining  these  pieces  together,  the  third  order  Taylor  series 
expansion  is  obtained 
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Rn  -  ^ao  *  ^ - ; -  ■•■  —  -  ^  ^6^  [cos{2nP) 

^no  i?no 

i?ii?2sin2{273P)  1  i2ii?2  53  r  -3^i;?2Cos(2np)sin(2^iP)  (31) 

d/ 2  -^  60/  771  " 


-  sin 


(2„p)*  3i;A=sinM2nP)  ^ 


(2)  Phase  Approximation.  R\*  and  R'^'  play 
significant  roles  in  the  phase.  Placing  these  terms  in  the 
phase  in  the  third  order  Taylor  series  expansion  and 
rearranging  the  terms  gives 


exp(jVt)expi  (-jic)  [ {rL)  ±  5  cos(2np) 


Rj^R^sin^  (2n^) 


'/lO 


■^1^2  g3  -3i?ii?2<^®®  (2nP)  sin  (2nP) 

SRL  Rr  ^ 


-  sin (2np)  + 


SR^^JJj^sin^  (2nP)  ""i  2 


r'  * 

^no 


->5 


(32) 


Since  e***  does  not  contain  the  index  n,  it  can  be  set  aside  for 
the  moment  and  the  second  term  will  be  examined  for  possible 
reductions.  In  order  for  a  term  to  be  considered  negligible, 
a  phase  error  of  less  than  7r/4  is  assumed.  Certain 
assumptions  also  need  to  be  made  in  order  to  represent  a 
typical  ocean  environment.  First,  it  is  assumed  that 


iO.l 


which  is  equivalent  to 


(33) 
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(34) 


sin0p:SO  ,  4 

Where  c,  is  the  speed  of  sound  in  the  water,  Cj  is  the  speed  of 
sound  in  the  bottom,  and  6^  is  the  critical  angle. 

Second, 

piO.l 


Finally, 


(36) 


Whose  implementation  will  be  seen  in  the  following  analysis. 

Grouping  the  terms  inside  the  brac)cets  from  Eq.  32  in 
terms  of  powers  in  R,  and  Rj  gives 
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(37) 


,  ^1^2  r  1 

^—TT  t|sin3(272P)63] 

^no 

[^sin2(272p)62-|cos{222P)sin(2i2p)63]  (38) 

^no 

R  R  1 

[sin(2/3p)  5  +  -icos  (2np)  52- Agin (2np)  6^]  (39) 

Rno  2  ® 

kRL  <40) 

Each  group  needs  to  be  compared  to  7r/4  to  determine  if  it  can 
be  considered  negligible.  Beginning  with  expression  37  the 
following  inequality  is  set  up 

[^sin3(2hP)62]<it/4  (41) 

^no 

Rearranging  gives 

k(~)  (-^)  i?2  [-|sinM2nP)  6^]<7i/4  (42) 

p'  p'  z 

^170  ^no 

Making  the  definitions  S=Ri/R'^  and  sinn=R2/R'no  multiplying 
through  by  X/X  gives 

kX  5^sin2Q-^  [isin^  (2nP)  6^]  <n/4  (43) 

where  kX=7r/ (2sin^5tan/3)  for  a  fast  bottom. 

To  ensure  that  the  left  hand  side  of  the  equation  is  always 
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«  77/4  the  maximum  of  all  powers  of  sine  and  cosine  terms  will 
be  taken  to  be  unity,  with  the  exception  of  sin^,  which,  as 
stated  in  the  above  assumptions,  already  has  a  maximum  value 
of  0.4.  Substituting  these  maximum  values  into  Eq.  43  and 
rearranging  gives 


^  53  ^  ^  JL 

.4tanp  X  53 


(44) 


Since  /8  is  small,  tan  =  /3  while  8,  the  receiver  angle,  falls 
between  0  and  jS.  The  right  hand  side  is  minimized  by  setting 
^=8.  Making  this  substitution  into  Eq.  44  and  allowing  to 
take  on  a  minimum  value  of  0.1,  reduces  Eq.  44  to 


i<40 


(45) 


Finally,  substituting  S  =  1/10  according  to  Eq.  36  gives 


^<40000  (46) 

If  Rj/X  for  some  reason  became  larger  than  40000,  the  term 
would  need  to  be  added  back  into  the  equation  to  keep  the 
phase  error  less  than  7r/4.  Note  that  j8  >  0.1  is  eqivalent  to 
setting  a  lower  limit  on  the  wedge  angle  of  about  6  degrees. 
Expression  38  can  be  written 
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(47) 


t4sin2(2np)fi' 

p/  nf  ^ 

i^riO 


-icos(2nP)sin(2np)5^1  <it/4 
2 

„a^in,  tne  sa.e  substitutions  as  in  the  ptsvious  aquation 
yields 


JH _ 52-^  Csin2  (2np)  52 


"aslne^tanp  X 

-  cos  (2jaP)  sin (2np)S^l  <11/2 

ones  again  substituting  in  the  saxisising  values  gives 


(48) 


52 -^<4 
X 


(49) 


substituting  S  *=  l/lO  gives 


-?a<400 

X 


(50) 


which  allows  this  expression  to  be  neglected. 


Expression  39  gives 

Ar,  lsln(2nP) « 4c°s  (2"P'  SMef 


(51) 


k — ^*^2 


Making  the  same 


substitutions  as 


in  the  previous  equation 


yields 
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(3)  Amplitude  Approximation.  It  can  be  shown 
that,  in  the  far  field,  the  amplitude,  R'„*  and  R'„'  can  be 
replaced  by  R'^  with  a  relatively  small  margin  of  error. 
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and  sina 


b.  Derivation  of  sina'„'^ 

To  take  the  Taylor  series  expansion  of  sina',/  and 
sina ,  an  expression  for  each  term  is  needed.  The  geometry 
for  this  is  shown  in  Fig.  li  for  an  upper  doublet. 


Geometrical  Representation  of  *  and  R„’* 
Figure  11 

From  the  law  of  sines 


sino^* 

Rearranging  this  gives 

R'*3ino'„*=sinr)^I?2^*yJ 

Substitution  Vl-cos^  17  for  sin  tj  gives 

i?'"sina^,"  -  {R^^*Y^)  cos'^r\ 


(56) 


(57) 


(58) 


To  get  T)  in  terms  of  5  a  substitution  for  cos^  ij  is  needed. 
To  do  this  the  following  equality  is  setup; 
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i?i  ^rI  +  Y0-2 Ri  +  Y^cosr\ 

Canceling  similar  terms  and  solving  for  cos^  t/  gives 

2  l?|cos2  (2np+6) 

cos^n  =  — = - !- 1  (60) 

i?2+1^0 

Substituting  this  into  Eg.  58  and  replacing  cos^(2n|8  +  S) 
with  l-sin^(2n/3  +  S)  gives: 

i?n^sino^*=^y/+i?2^sin^  (2nP+6) 


Bringing  R'„*  over  to  the  right  side  gives  the  final 
expression  for  sina'^*. 


And  finally, 
is 


sina„  = 


♦  _  y/y^^+i?2^sin^ (2np+6) 


R 


/* 


(62) 


by  the  same  reasoning  the  expression  for  sina 


sino 


3 in^  (2np  -6 ) 


Rr 


(63) 


When  the  receiver  is  moved  vertically  to  the  surface  Eq.  62 
and  Eq.  63  become 


sino. 


_  ylYj^*R2^sin^(2n^) 

rL 


Again,  by  setting  Yo=0  for  a  geometrical  check 


(64) 
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(65) 


s  ina„*  = 


Rr 


^^2sin(2iip-5) 


(66) 


which  agrees  with  the  upslope/downslope  case. 

(1)  The  Taylor  series  expansion  for  sina  and 
sina\'.  The  first  order  Taylor  series  expansion  taken  around 
Yo=0  gives 


,  .  ,  -  dsinaj . 

sino>sino„^'+5 - l».o 


(67) 


Because  these  terms  only  appear  in  the  amplitude,  the 
approximation  sina »  sina'no  is  acceptably  accurate. 

c.  Total  Pressure  Including  Approximations 

Putting  all  of  these  equations  together  yields  the 
pressure  for  the  upper  doublet, 


Pn  [sin^k^dsinoao] 

Rno  2 

-jJci  (sin(2nP)#*-lco8(2nP)ft*)  J 

^  rL  2 


(68) 


That  for  the  lower  doublet  is  the  same  except  for  a  sign, 
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4.  Validation 

To  validate  the  analysis  of  doublet  cross-slope 
propagation,  a  computer  program  was  developed.  See  Appendix 
"A”.  Different  inputs  of  Y^,  R,,  and  Rj  were  used  along  with 
a  wedge  angle  and  source  depth  of  7.5  and  3.75  degrees 
respectively.  These  inputs  were  also  put  into  the  program 
URTEXT  and  a  comparative  table  was  constructed.  See  Table  1. 
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table  1 


1  COMPARISON 

OF  DOUBLEr  ANALYSI5 

;  TO  IMAGE  I’HEORY  1 

p  =  7.5  Y  = 

3.75 

#  of  depth  points 

=  6 

R1  = 

1  R2  = 

15  YO 

=  0 

R1  = 

1  R2  = 

12  Y( 

)=  8 

1  DOUBLEr 

!MAGE 

DOUBLET 

IMAGE  1 

AMP 

P!!S 

AMP 

PHS 

AMP 

PHS 

AMP 

PHS 

.19230 

-1.189 

.19020 

-1.174 

.20760 

.97186 

.14901 

.99135 

.37155 

-1.192 

.36747 

-1.175 

.40214 

.97103 

.28788 

.99132 

.52564 

-1.195 

.51973 

-1.177 

.57139 

.96959 

.40713 

.99133 

.64414 

-1.199 

.63659 

-1.179 

.70471 

.96741 

.49862 

.99135 

.71913 

-1.203 

.71006 

-1.180 

.79371 

.96432 

.55614 

.99135 

.74575 

-1.207 

.73519 

-1.180 

.83276 

.96002 

.57590 

.99138 

AAMP  = 

=  .0066 

APHS  =  .01898 

AAMP 

=  .1729 

APHS  = 

.0232 

R1  = 

1  R2  = 

14  YO 

=  3 

R1  = 

=  1  R2 

=  8  YO 

=  9 

1  DOUBLE!' 

IMAGE 

DOUBLEr  1 

IMAGE  1 

AMP 

PUS 

AMP 

PHS 

AMP 

PHS 

AMP 

PHS 

.16437 

.77749 

.15456 

.80884 

.18308 

-.2626 

.12944 

-.2403 

.31909 

.77590 

.29861 

.80795 

.35338 

-.2536 

.25005 

-.2402 

.45507 

.77353 

.42229 

.80677 

.49899 

-.2554 

.35361 

-.2402 

.56437 

.77059 

.51721 

.80558 

.60956 

-.2580 

.43309 

-.2403 

.64061 

.76720 

.57687 

.80473 

.67701 

-.2617 

.48304 

-.2403 

.67940 

.76316 

.59736 

.80443 

.69602 

-.2671 

.50009 

-.2403 

AAMP 

=  .0423 

APHS  =  . 

0351 

AAMP 

=  .1447 

APHS  = 

.0178 

R1  = 

1  R2  = 

13  YO 

=  6 

R1  = 

1  R2  = 

o 

=  10 

1  DOUBLE!' 

IMAGE 

DOUBLET 

IMAGE  1 

AMP 

PHS 

AMP 

PHS 

AMP 

PHS 

AMP 

PHS 

.19161 

.77072 

.14813 

.74146 

.00110 

-.3287 

.00016 

-.0074 

.36918 

.77044 

.28616 

.74157 

.00216 

-.3376 

.00030 

-.0376 

.51968 

.76982 

.40469 

.74174 

.00320 

-.3477 

.00043 

-.0112 

.63207 

.76894 

.49564 

.74189 

.00415 

-.3607 

.00052 

-.0134 

.69820 

.76776 

.55281 

.74199 

.00500 

-.3813 

.00055 

-.0195 

.71336 

.76620 

.57226 

.74203 

.00573 

-.4041 

.00055 

-.0251 

AAMP 

=  .11073 

APHS  = 

.0272 

AAMP 

=  .0031 

APHS  = 

.3421 
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It  is  shown  that  the  doublet  cross-slope  propagation  results 
are  in  good  agreement  with  the  URTEXT  results.  It  should  be 
noted  though,  that  for  the  results  to  compare  favorably,  the 
wedge  angle  that  is  chosen  must  ensure  that  the  number  of 
wedges  are  an  even  multiple  of  four. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


Image  theory  is  an  excellent  method  for  predicting  the 
transmission  of  sound  in  shallow  water.  It  is  an  important 
area  of  research  because  it  is  directly  applicable  to  shallow 
water  Anti-submarine  Warfare.  Other  methods  that  are  being 
examined  include  the  Parabolic  Equation  approximation,  ray 
tracing  methods,  and  the  Adiabatic  Normal  Mode  theory.  Image 
theory  shows  the  most  promise  for  predicting  the  acoustical 
pressure  fields  in  a  wedge-shaped  ocean  because  it  can  predict 
pressure  fields  for  cross-slope  and  for  all  frequencies. 
Image  theory  can  also  deal  with  the  transitional  area  at  cut 
off,  a  region  that  is  difficult  for  normal  mode  theory. 

Previous  theses  have  laid  the  groundwork  for  the  image 
theory  model.  This  thesis  has  extended  the  doublet 
approximation  to  cross-slope  propagation.  The  base 
approximation  used  ,  R,'  «  limits  the  extension  to  the 
far-field.  Another  limiting  restriction  on  this  thesis 
basically  states  that  for  a  phase  error  of  7r/4  or  less,  the 
scaled  receiver  distance  R2/X  must  be  less  tha,!  400. 
Geometrically  this  says  that  when  the  receiver  gets  too  far 
away  the  phase  errors  combine  and  the  pressure  equation 
developed  in  Eq.  72  is  no  longer  valid.  One  further 
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restriction  is  that  the  wedge  angle  must  be  chosen  such  that 
the  number  of  wedges  produced  is  an  even  multiple  of  four. 
Because  of  this,  when  the  dipoles  are  grouped  into  pairs  of 
doublets,  one  dipole,  180  degrees  from  the  neutral  doublet, 
will  be  left  ungrouped.  This  will  only  be  a  significant 
factor  when  the  grouped  dipoles  cancel  completely  leaving  only 
the  dipole  to  contribute  to  the  pressure.  The  effects  of  this 
dipole  can  also  be  seen  in  Table  1:  The  results  vary  when  Rj 
gets  small  and  gets  large.  Although  this  dipole  can  have 
noticeable  effects  on  the  pressure  in  the  rigid  bottom  case, 
when  reflection  coefficients  are  incorporated,  the  effects 
become  negligible. 

The  next  step  in  this  research  would  be  to  incorporate  the 
reflection  coefficients  for  the  upper  and  lower  doublet  into 
Eq.  72.  The  reflection  coefficients  though,  because  of  their 
large  effect  on  the  pressure,  need  to  be  calculated 
individually  for  each  image.  Nassopoulos  incorporated  the 
reflection  coefficients  into  the  upslope/downslope  case,  but 
an  extension  is  still  needed  for  the  cross-slope  case. 
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APPENDIX  ''A” 


1''  lUM  t-T"''!:  L>>:fr>ir  ,  I'lrol  r  ^^^f■Rf'XIHATION 

20  PV.M  cr-?"-?  slono.  witl^  rigid  bottom.  Including  third  order  approx 

"in  pt  ,  l-l  '  '  '2'  N 

iiO  fflf't'T  "  =inn]  osf  )  nf  ,.jodoe  .  source;  B.G  =*,B,G 
Z-0  INPUT  "Ml  =  B2(deU.a  P) 

60  MClfj  INH  ?60  Br  .oool  )  2  )/4  )+ .0001  ) 

70  INPUT  "source,  receiver,  apex  distances:  R1,R2,Y0  a',Rl.R2,Y0 
eo  INPi.T  ■‘‘■l.'rS 
OO  rt'ieo/pi:  B=B.-'T6;  G=G/T6 

IOC  STORCO  ?  -  ?QP' 0B'^(  C*C  ■  1  )  ):  T4;;Pl/(  2»STARCOS«-TAN(  B  )  ) 

110  LPRINT  U'rING  •'THPSIS  APPROX  angles:  wedge*  «#.*  source*  m» .  * "  ;B*T6  .G*T6 
120  l.pRiru  UOING  "distances:  sour  ce*!*#  .#11  receiver*##.##  shore*## -#*  ;R1  .R2  .YO 

1  .’0  LPRINT  USING  "  cl/c2*  #.««  "iC 

lilC  LPPIMT  "SING  "k:''-pi/(r  tenB  sgr  abs(  1-c' 2  ))*«#«.####":  T4 
jr-o  LPRINT  USING  "No.  of  Dipole-Pairs  =#«.#«'SNPTS 
If  O  Li  PIN  I  ‘  DFG  APP  PHASE  " 

ISS  niST^2rRl*SlN(  G  ) 

170  Y2^YO*vo:  P22*R2'P.? 

ISO  i:2*C'*C:  D2*Y0»Y0*Rt*Rl  iR2*R2 
IOC  I’ U'2 

200  roR  M-il  ro  inn;  U  IMBcmi:  IF  (B-D)<-.0001  THEN  3S0 
210  R0*S0R(  P2 -RS*C  0S(  P  )  ):  GAM0-'R1*R2/R0 
220  I:”!  2/RO  irSINl  M ♦GAPOrG'-D  ) 

220  Pl*r*SlN(  T4*PO  ):  P2*r«COS( T4*R0  ) 

2-10  P2*0:  P-l-iO 

2S0  FOR  N«1  10  NPTS 

2SC  RN*SQR(  1'2  RS*L0S(  2  UirB  )  ) 

262  RNl*(T4»Rl»R2tSlN(2rN*8))/RN 
264  RN2«( T4»Rl*R2rC06(  2*N*B))/RN 
270  jNS»(S0R(  Y2*R22«(  SIN(  2*N*B  )  )''2  ))/RN 
20C  X"i(('l  N  ;♦(  2/RN  )rGIN(  .  S*T4+DIST*SNS  ) 

2B2  THDPOS*  -T4rRN  -  D*RN1  -  .5*0"2-»RN2 

202  rHDNEG-  -T4rRN  +  D*PN1  -  .Srp'2*RN2 

2B4  THOR  *  C0S( THDPOS  )-C0S(  THDNEG) 

296  THPI  S1N(  thppOS  ) -S1N(  THONEG  ) 

2B7  REALl  =  -X*THDI 

2«S  IMAGl  -  X-»THPR 

200  PO-’REALl  »  P? 

210  pa-TMAGl  •  P4 
220  NEXT  N 

2  20  PR-P11-P2:  PIM-^2^P4 
240  AMP-SOR(PR*PR»PrM*PIM  ) 

?SC  PHS*A1N(  PIM/PR ) 

260  LPRINT  USING  "  ##.#«  ##.#«#««  #«.#####  D*T6,AMP.PHS 

3,’0  NEXT  M 
290  END 
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